Multiband model for tunneling in MgB2 junctions 
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A theoretical model for quasiparticle and Josephson tunneling in multiband superconductors is developed and 
applied to MgB2 -based junctions. The gap functions in different bands in MgB2 are obtained from an extended 
Eliashberg formalism, using the results of band structure calculations. The temperature and angle dependencies 
of MgB2 tunneling spectra and the Josephson critical current are calculated. The conditions for observing one 
or two gaps are given. We argue that the model may help to settle the current debate concerning two-band 
superconductivity in MgB2. 
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Soon after the discovery of superconductivity in MgB2,u 
first principle calculations were performed to determine the 
electronic structure of this material. It was found that the 
Fermi surface consists of two three-dimensional sheets, from 
the TT bonding and antibonding bands, and twcupsarly cylindri- 
cal sheets from the two-dimensional a bands.fflffl The multi- 
band picture has given rise to tbe^oncept that two supercon- 
ducting energy gaps can coexistfltl in MgB2. 

Two-band superconductivity is. a phenomenon that has been 



E£cent experimental STM 
lil high-resplution photo- 



observed in Nb doped SrTiOa. 
and point-contact specjtipscopy^ 

emission spectmscopyjlj Raman spectroscopy ,0 specific heat 
measurementstj and muon-spin-relaxation studies of the mag- 
netic penetration depthO support the concept of a double gap 
in MgB2 (see Ref.|l^ for a review of experiments). However, 

itpoint-contact data 
Moreover, some 



there is an ambiguity in the interpretation-a3 
concerning the existenc£-pf two gaps.&t^ 



tunneling measurementsU show only one gap with a magni- 
tude smaller than the BCS value of A = 1.76 ksTc- 

In order to resolve this discrepancy, we address the ques- 
tion, how multiband superconductivity will manifest itself in 
tunneling. We present the theoretical model for quasiparticle 
and Josephson tunneling in MgB2-based junctions. Using the 
results of band-structure calculations, we apply an extended 
Eliashberg formalism to obtain the gap functions in different 
bands, taking strong coupling effects into account. Tunneling 
from a normal metal (N) into MgB2 is considered, in an ex- 
tended Blonder- Tinkham-Klapwijk (BTK) model .Ei] The tem- 
perature dependencies and absolute values of the IcRn prod- 
uct (Ic is the critical cuiTent and Rn is the normal state re- 
sistance) are calculated in MgB2-based SIS tunnel junctions, 
where S denotes a superconductor and I an insulator Tun- 
neling in the direction of the a-b plane, in the c-axis direc- 
tion and under arbitrary angle is considered. Furthermore, the 
Josephson supercurrent between a single-gap superconductor 
and MgB2 is calculated. 

According to the labeling of Liu et al.p the four Fermi sur- 
face sheets in MgB2 are grouped into quasi-two-dimensional 
a bands and three-dimensional tt bands. Hence, normal and 
superconducting properties of MgB2 can be described by an 



effective two-band model. Within this model, Liu et alE es- 
timated the coupling constants and energy gap ratio in the 
weak coupling regime. More recently, the band decompo- 
sition of the superconducting and transport Eliashberg func- 
tions afjFij (lli) (where i and j denote cr or tt bands), which 
describe the electron-phonon coupling in MgB2 as function of 
the frequency uj, was provided in Ref.[l^. This allows to per- 
form a strong coupling calculation of the superconducting en- 
ergy gap functions (w„) in different bands. The functions 
Ai {ojn) in turn determine the Josephson critical current in a 
tunnel junction between multiband superconductors, which is 
giveruhy a straightforward generalization of thei-well-known 
resultEj to the case of several conducting bandsEll as well as 
strong coupling. The critical current component for tunneling 
from band i into j is given by 
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(1) 



where C and TZ denote left and right superconductors respec- 
tively, R^j^ = min{i?^j?^, R^lj} is the normal-state conduc- 
tance of a junction for the bands (i, j) which is given by the 
integral over the Fermi surface SciCRj) 
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where A is the junction area, w„ is the projection of the Fermi 
velocity vf on the direction normal to the junction plane, and 
Dij is the probability for a quasiparticle to tunnel from band i 
in C into band j in TZ. The total critical current is the sum of 
the components Ic = Tlij ^v- 

The gap functions A^ can-be calculated with an exten- 
sion of the Eliashberg formalismEl to two bands 



+ A2 {to, 



(3) 
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A2 {lj„,) 



(4) 



^af^Fii{uj)du!. The cor- 
= 0.042, 
are 



where Ay 

(wn) are the Migdal renormalization functions and uj„ 
7rT(2n + 1). These equations are solved numerically with 
the electron-phonon a^jFij [oj) functions from Ref.|T^. The 
cutoff frequency lOc is taken equal to 10 times the maximum 
phonon frequency. The functions /I*^ represent a matrix of 
the Coulomb pseudopotentials defined at lOc, calculated up 
to a common prefactor that is used as an adjustable param- 
eter to get Tc = 39.4 K. The matrix ji* at the frequency uj\a 
(relevant for the McMillan expression for Tc in the isotropic 
case) is given by /i* = [1 + Ji* hi{ijjc/ ^jJ\n)\^^V* , where a;i„ 
follows from = h\{Lj/L0\a)L0^^a1jFij 
responding matrix elements are /i*^. = 0.13, /i*^ 
/i*^ 0.03, iil^ = 0.11 and Ay = w„) from Refjl9 
\aa = 1-017, = 0.213, = 0.155, A^^ = 0.448. Due 
to the interband coupling terms in Eqs. ^ both gaps close 
at the same Tc- The resulting temperature dependencies of the 
energy gaps, (T) , are plotted in the inset of Fig. |l] and it 
is found that (T = 0) = 7.09 meV and A^ {T ^Q) = 2.70 
meV, with the 2A/Tc ratios being equal to 4.18 and 1.59, re- 
spectively. For comparison, also the BCS curve is shown for 
Tc = 39.4 K. The BCS value for the gap that corresponds to Tc 
= 39 A K is 6.0 meV at K. It can be seen that the temperature 
dependencies are qualitatively different from the BCS temper- 
ature dependence. The ratio of the gaps Acr/Ajr increases for 
increasing temperatures, as was experimentally observed for 
example in Ref.^. 

The influence of impurities can be incorporated into the 
model. Intraband scattering does not change the two gaps 
(Anderson's theorem), while the interband scattering can be 

included by terms -y.jAj/Juj'^ + A|, ^ijUn/ ^uj^ + A^ in 

the Eliashberg equations (|^, ^ respectively. The smallness 
of 7y compared to ttTc indicates that the double-gap feature 
should experimentally be observable, also in thin films, even 
for a certain amount of impurity scattering. A large amount 
of impurity scattering (7^ exceeding the maximum phonon 
frequency) will cause the gaps to converge to the same value. 
From Eqs. (^ ^ and including the scattering terms an asymp- 
totic value of Act = Ajr = 4. 1 meV and Tc ==25.4 K is found, 
giving a 2/S./kBTc ratio of 3.7. 

In order to obtain the normal state resistance, we have to 
evaluate the effective junction transparency components Dij 
In the case of a specular barrier, U {x) — Uq5{x — xq) 
given by 



D^j IS 



D, 



J {VnXi + Vn.TZj) +UQ/h 



(5) 



It follows from Eqs. (^ g) and as first pointed out in Ref.|22|, 
that the normal state conductance R^^^ in the large Uq limit is 
proportional to the Fermi-surface average (^Nv^) ^. The latter 
is proportional to the contribution of the electrons in band i to 
the squared plasma frequency (w*) . This essentially simpli- 
fies the task of summing up the interband currents since the 



TABLE I: Calculated plasma frequencies, average Fermi velocities 
and gap values for the o and vr bands. 





(eV) 


Lol (eV) 




u^(m/s) 


A (meV) 


(7 band 


4.14 


0.68 


4.40 10'' 


0.72 lO'' 


7.09 


TT band 


5.89 


6.85 


5.35 lO'^ 


6.23 lO'^ 


2.70 



partial plasma fremiericies are available from the band struc- 
ture calculations fflcra The normal state junction conductance 

is thus proportional to ((.^p)^ ('^n)7^J■' where {vn)-jzj is the 
average Fermi velocity projection in the corresponding band 
(see Table ||). In order to sum up the contributions of differ- 
ent bands, we restrict ourselves to the weighing factors (w*) , 
neglecting the difference in {vn)-jij- This is a reasonable ap- 
proximation since the difference between vp in the a and tt 
bands in the a-b plane is rather small, while for c-axis tunnel- 
ing only the tt band contributes, as will be shown later, so that 
the problem of summation does not appear in this case. 

SIN tunneling. The conductance in a MgB2-I-N tunnel 
junction is the sum of the contributions of twp-jbands. Each of 
the conductances is given by the BTK model,lUl where the cor- 
responding normal state conductances are proportional 
to the minimum of the square of the plasma frequencies at the 
N and MgB2 sides. Since the plasma frequency in a typical 
normal metal (e.g. Au, Ag) is larger than the plasma frequen- 
cies in MgB2, the conductances are limited by the electrons 
on the MgB2 side 



(6) 



Finally, the normalized conductance of an N-I-MgB2 contact 
is given by 



aiV) 



\dV ) 



NIS 



' dl \ 

•.dVI NIN 



<j\2 



(7) 



Here, the dimensionless conductances (Ja.-K {V) are provided 
by the BTK model, with the calculated values for the gaps and 
plasma frequencies, as shown in Table ^ 

In the conductance versus voltage plot (Fig. ^ for tunneling 
in the a-b direction, two peaks are clearly jtfisible, in qualita- 
tive agreement with the experimental dataBtUiil The ratio of 
peak magnitudes is not only determined by the ratio of the 
plasma frequencies, but also by thermal rounding and by the 
barrier strength Zbtk — Uo/hvp (where Wi? is taken constant 
for the different bands for the same reason as was given in the 
determination of Rij). In particular, the peak at the smaller 
gap dominates in the small Zbtk regime (point-contact), while 
the second peak dominates at large values of Zbtk (tunnel- 
ing), as may be seen in Fig. |l| at 4.2 K. 

Due to the smallness of oj^ in the c-direction, it can be seen 
from Eq. (^ that the conductance in the c-direction is only 
determined by the tt band. In this case, no double-peak struc- 
ture is expected in the conductance spectrum. This explains, 
together with the dependence on Zm^K, why in some experi- 
ments only.pne.peak was observedliJ or why the second peak 
was weak. 



one. 
lojii 
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V(mV) T{K) 



FIG. 1: The normalized conductance of MgB2-I-N junctions as 
function of voltage at 4.2 K, with Z — Uo/hvp- The supercon- 
ducting gaps for the different bands are shown in the inset. 



FIG. 2: IcRn temperature dependence for different tunneling com- 
ponents. The resulting IcRn for tunneling in the direction of the 
a-b plane and c-axis direction are indicated, for MgB2-I-MgB2 and 
Nb-I-MgB2 junctions (see text). 



Note, that the assumption of the ratio of the normal state 
conductivities being equal to the ratio of the square of the 
plasma frequencies holds when the interface is a (5-function 
shaped tunnel barrier, with large Z^tk- This means that for 
small Zbtk, the results should be considered as a qualitative 
indication only. In the latter case, as well as for other types 
of barriers, a numerical integration of Eqs. (||, ||) must be per- 
formed. 

SIS Josephson tunneling. We consider Josephson tunnel- 
ing between two MgB2 superconductors. With the values for 
the plasma frequencies, ujp < ujp, this gives Ra-K = R-wa — 
ma-x {Raa , Rtttt) = Raa and R^T^/R^^ = R^^/R^^ — 
{ujp/uipY > 1- The total conductance is given by i?^^ = 

ED-l 
ij ^ij ■ 

For tunneling in the a-b plane (as can be realized for ex- 
ample in an edge configuration), with Ra-^r — Rttct and 



/ttct, the total IcRn product becomes 



IcRn — 



IaaR(7(T ~^ '^la-K Ran ^ Itttt Rtttt (^p / ) 



(8) 



The results of numerical calculations are presented in Fig. |[ 
Due to strong-coupling and interband coupling effects, the 
temperature dependencies of lijRij differ from the well- 
known Ambegaokar-Baratoff result for an SIS junction be- 
tween isotropic superconductors, most clearly demonstrated 
by the positive curvature of the I^^R^^^ contribution. The 
IcRn value at T = 4.2 K is 5.9 mV. 

For tunneling along the c-axis, the only contribution to the 
IcRn product comes from IttttRtttt, because of the negligible 
value for ojp in the c-axis direction. This gives IcRn = 4.0 
mV at T = 4.2 K. 

The plasma frequency in a certain direction, under an an- 
gle ip with the a-b plane, can be determined from the ellip- 
soid equation {ujj,)^ = {ujI^^ + {oj\,,^f , and (w^ ^_ 



(cjp^/Wp^^j)^ = 1, where and cj^^ form the decom- 
position of Because of the negligible value of w^^' 
it is evident that is negligible for nonzero values of 

Lp = arctan((x'p ^/'^p x)- This implies that tunneling under 
a nonzero angle with the a-b plane gives the same result as 
tunneling in the c-axis direction, namely IcRn = 4.0 mV at 
T = 4.2 K. For angles approaching zero (of the order of 0.6°), 
IcRn rapidly increases towards the maximal value for tunnel- 
ing from a-b plane to a-b plane, namely IcRn = 5.9 mV at T 
= 4.2 K. For a large amount of impurity scattering the lijRij 
values converge to the same value. It follows in that case from 
Eq. (^, with the plasma frequencies from Table |[ that IcRn 
becomes almost isotropic. 

Finally, tunneling from MgB2 into a superconductor S' 
with a single gap will be considered (we take Nb as an ex- 
ample). The resulting liS'RiS' temperature dependencies are 
calculated numerically, using 1.4 mV for the energy gap in 
Nb. The ratio of resistances is determined from Eq. (S). Since 
typical values of plasma frequencies in other superconductors 
are bigger than in MgB2 (e.g. 9.47 eV for Nb, 12.29 eV for 
Al and 14.93 eV for Pb, see Ref.|24|), the following expression 
is obtained 



IcRn 



hs'RaS'+hs'R^S' (Wp/Wp) 



1 



(9) 



when tunneUng occurs into the a-b plane of the MgB2. In the 
case of c-axis tunneling, only the /^rS' RirS' contribution re- 
mains. The results for tunneling from Nb to MgB2 are also 
indicated in Fig. ^. Other superconductors give qualitatively 
similar results. The only scaling parameter is the critical tem- 
perature of the superconducting counter-electrode. 

Our results for Josephson tunneling provide an upper bound 
for IcRn products, being 5.9 mV and 4.0 mV for tunneling 



4 



into the a-b plane and c direction respectively. There have 
already been several observations of Josephson currents in 
MgB2 junctions^ with IcRn values that are much lower than 
our predictions. This can be due to extrinsic reasons such as 
a degradation of the Tc of surface layers in the vicinity of the 
barrier, the barrier nature and barrier quality. From our model, 
however, it follows that polycrystallinity does not reduce the 
Josephson coupling very much, as indicated by the calculated 
value of IcRn of 4.0 mV for c-axis transport, neither does 
strong impurity scattering because of the relatively large aver- 
age gap of 4. 1 meV in this case. 

In conclusion, Josephson tunneling in MgB2-based junc- 
tions is discussed theoretically in the framework of a two- 
band model. The gap functions in different electronic bands 
are calculated using the Eliashberg formalism together with 
band structure information. This provides a basis to interpret 
electronic transport in MgB2. We have shown the possibil- 



ity to observe either one or two gaps in point-contact spectra 
of MgB2, depending on the tunneling direction, barrier type 
and amount of impurities. The results are also relevant for the 
electronic application of MgB2 since they provide the limit 
for the Josephson coupling strength in MgB2 based junctions. 
For MgB2 in the clean limit we have shown that IcRn values 
as high as 5.9 mV can be expected for MgB2 tunnel junctions 
if tunneling occurs in the direction of the a-b plane. In other 
cases the limiting IcRn values will not exceed 4.0 mV. Our 
predictions for the gap and IcRn anisotropy and for the Ic 
vs T dependence in MgB2-based junctions can be verified ex- 
perimentally and thus may help to settle the current debate on 
two-band superconductivity in MgB2. 
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